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Summary. We introduce the snowblower problem (SBP), a new optimization problem that is 
closely related to milling problems and to some material-handling problems. The objective in 
the SBP is to compute a short tour for the snowblower to follow to remove all the snow from 
a domain (driveway, sidewalk, etc.). When a snowblower passes over each region along the 
tour, it displaces snow into a nearby region. The constraint is that if the snow is piled too high, 
then the snowblower cannot clear the pile. 

We give an algorithmic study of the SBP. We show that in general, the problem is NP- 
complete, and we present polynomial-time approximation algorithms for removing snow un- 
der various assumptions about the operation of the snowblower. Most commercially-available 
snowblowers allow the user to control the direction in which the snow is thrown. We differ- 
entiate between the cases in which the snow can be thrown in any direction, in any direction 
except backwards, and only to the right. For all cases, we give constant-factor approximation 
algorithms; the constants increase as the throw direction becomes more restricted. 

Our results are also applicable to robotic vacuuming (or lawnmowing) with bounded ca- 
pacity dust bin and to some versions of material-handling problems, in which the goal is to 
rearrange cartons on the floor of a warehouse. 



1 Introduction 

During a recent major snowstorm in the northeastern USA, one of the authors 
used a snowblower to clear an expansive driveway. A snowblower is a "material 
shifting machine", lifting material from one place, and depositing it nearby. There 
is a skill in making sure that the deposited piles of snow do not grow higher than 
the maximum depth capacity of the snowblower. This experience crystallized into an 
algorithmic question we have called the Snowblower Problem (SBP): How does one 
optimally use a snowblower to clear a given polygonal region? 
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(NAG2-1620), NSF (CCR-0098172, ACI-0328930, CCF-0431030), and Metron Aviation, 
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The SBP shows up also in other contexts: Consider a mobile robot that is 
equipped with a device that allows it to pick up a carton and then place the car- 
ton down again, possibly on a stack of cartons, in a location just next to it; with each 
such operation, the robot shifts a unit of "material". The SBP models the problem 
in which the robot is to move a set of boxes to a specified destination in the most 
efficient manner, subject to the constraint that it cannot stack boxes higher than a 
capacity bound. 

In a third motivating application, consider a robotic lawnmower or vacuum 
cleaner that has a catch basin for the clippings, leaves, dust, or other debris. The 
goal is to remove the debris from a region, with the constraint that the catch basin 
must be emptied (e.g., in the compost pile) whenever it gets full. 

The SBP is related to other problems on milling, vehicle routing, and traveling 
salesman tours, but there are two important new features: (a) material must be moved 
(snow must be thrown), and (b) material may not pile up too high. 

While the SBP arises naturally in these other application domains, for the rest of 
the paper, we use the terminology of snow removal. 

The objective of the SBP is to find the shortest snowblower tour that clears a 
domain P, assumed to be initially covered with snow at uniform depth 1 . An impor- 
tant parameter of the problem is the maximum snow depth D > 1 through which 
the snowblower can move. At all times no point of P should have snow of greater 
depth than D. The snow is to be moved to points outside of P. We assume that each 
point outside P is able to receive arbitrarily much snow (i.e., that the driveway is 
surrounded by a "cliff" over which we can toss as much snow as we want).^ 

Snowblowers offer the user the ability to control the direction in which the snow 
is thrown. Some throw directions are preferable over others; e.g., throwing the snow 
back into the user's face is undesirable. However, it can be cumbersome to change 
the throw direction too frequently during the course of clearing; also, some variants 
(e.g., lawnmowing with clippings) require that the throw direction is fixed. Thus, 
we consider three throw models. In the default model throwing the snow backwards 
is allowed. In the adjustable-throw model the snow can be thrown only to the left, 
right, or forward. In the fixed-throw model the snow is always thrown to the right. 
Even though it seems silly to allow the throw direction to be back into one's face, the 
default model is the starting point for the analysis of other models and is equivalent 
to the vacuum cleaner problem (discussed later). 

Results. In this paper we introduce the snowblower problem, model its variants, 
and give the first algorithmic results for its solution. We show that the SBP is NP- 
complete for multiply connected domains P. Our main results are constant-factor 
approximation algorithms for each of the three throw models, assuming D > 2; refer 
to Table 1 . The approximation ratio of our algorithms increases as the throw direction 
becomes more restricted. We give extensions for clearing polygons with holes, both 
where the holes are obstacles and cliffs. Then we discuss how to adapt our algorithms 



^ The "cliff" assumption accurately models the capacitated vacuum cleaner problem for 
which there is a (central) "dustpan vac" in the baseboard, where a robotic vacuum cleaner 
may empty its load [1] and applies also to urban snow removal using snow melters [2]. 
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for clearing nonrectilinear polygons and polygons with uneven initial distributions 
of snow. We conclude by giving a succinct representation of the snowblower tour, in 
which the tour specification is polynomial in the complexity of the input polygon. 
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Table 1: Approximation factors of our algorithms in the defauU model (left), adjustable-throw 
model (middle), and fixed-throw model (right). 



Related Work. The SBP is closely related to milling and lawn-mowing prob- 
lems, which have been studied extensively in the NC-machining and computational- 
geometry literatures; see e.g., [3,4, 10]. The SBP is also closely related to material- 
handling problems, in which the goal is to rearrange a set of objects (e.g., cartons) 
within a storage facility; see [8,9, 13]. The SBP may be considered as an interme- 
diate point between the TSP/lawnmowing/milling problems and material-handling 
problems. Indeed, for £> = oo, the SBP is that of optimal milling. Unlike most 
material-handling problems, the SBP formulation allows the material (snow) to pile 
up on a single pixel, and it is this compressibility of the material that distinguishes 
the SBP from previously studied material-handling problems. With TSP and related 
problems, a pixel only is visited a constant number of times, whereas with material- 
handling problems, pixels may have to be visited a number of times exponential in 
the input size. For this reason, material-handling problems are not even known to be 
in NP [8,9], in contrast with the SBP. Note that in material handling problems the 
objective is to minimize workload (distance traveled while loaded), while in the SBP 
(as in the milling/mowing problems) the objective is to minimize total travel distance 
(loaded or not). 

The SBP is also related to the earth-mover's distance (EMD), which is the min- 
imum amount of work needed to rearrange one distribution (of earth, snow, etc.) to 
another; see [7]. In the EMD literature, the question is explored mostly from an exis- 
tential point of view, rather than planning the actual process of rearrangement. In the 
SBP, we are interested in optimizing the length of the tour, and we do not necessarily 
know in advance the final distribution of the snow after it has been removed from P. 

Notation. The input is a polygonal domain, P. Since we are mainly concerned 
with proving constant factor approximation algorithms, it suffices to consider dis- 
tances measured according to the Li metric. We consider the snowblower to be an 
(axis-parallel) unit square that moves horizontally or vertically by unit steps. This 
justifies our assumption, in most of our discussion, that P is an integral-orthogonal 
simple polygon, which is comprised of a union of pixels - (closed) unit squares with 
disjoint interiors and integral coordinates. In Section 5 we remark how our methods 
extend to general (nonrectilinear) regions. 

We say that two pixels are adjacent or neighbors if they share a side; the degree 
of a pixel is the number of its neighbors. For a region R C P (subset of pixels). 
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let Gr denote the dual graph of R, having a vertex in the center of each pixel of R 
and edges between adjacent pixels. Sometimes when we speak of the region R, we 
implicitly mean the dual graph G^ . A pixel of degree less than four is a boundary 
pixel. For a boundary pixel, a side that is also on the boundary of P is called a 
boundary side. The set of boundary sides, dP, forms the boundary of P. We treat 
dP as a set of boundary sides, rather than just as a closed curve. We assume that 
the elements of dP are ordered as they are encountered when the boundary of P is 
traversed counterclockwise. 

An articulation vertex of a graph G is a vertex whose removal disconnects G. 
We assume that Gp has no articulation vertices. (Our algorithms can be adapted to 
regions having articulation vertices, at a possible increase in approximation ratio.) 

Algorithms Overview. Our algorithms proceed by clearing the polygon Voronoi- 
cell-by-Voronoi-cell, starting from the Voronoi cell of the garage g — the pixel on 
the boundary of P at which the snowblower tour starts and ends. The order of the 
boundary sides in dP provides a natural order in which to clear the cells. We observe 
that the Voronoi cell of each boundary side is a tree of one of two special types, 
which we call liyies and combs. We show how to clear the trees efficiently in each of 
the throw models. We prove that our algorithms give constant-factor approximations 
by charging the lengths of the tours produced by the algorithms to two lower bounds, 
described in the next section. 

2 Preliminaries 

Tie Breaking. Our rules for finding V{p) for a pixel p that is equidistant between 
two or more boundaries is based on the direction of the shortest path from p to V{p); 
horizontal edges are preferred to vertical, going down has higher priority than going 
up, going to the right — than going left. 

In details, let p be a pixel, equidistant from m boundary sides {ei...em} = fc 
dV. Let r{p), l{p), d{p), u{p) denote the pixels to the left, above, to the right and 
below p. First consider the case when p is a boundary pixel. If r{p) ^ P, then the 
right side of p is defined to be V{p); else if l{p) ^ P, the left side of p is V{p); else 
if d{p) ^ P, the bottom side of p is V{p); otherwise, the top side of p is V{p). If p 
is not a boundary pixel, we build the path fromp to V{p) as follows. If there exists a 
side e <E £ such that there exist a (shortest) path from p to e, going through r{p) then 
V{p) = V{r{p)); else if there exists a side e G £ such that there exist a (shortest) 
path from p to e, going through l{p) then V{p) = V{l{p)); else if there exists a side 
e E £ such that there exist a (shortest) path from p to e, going through d{p) then 
V{p) = V{d{p)); otherwise V{p) = V{u{p)). 

We first describe the rules for finding the closest boundary side for pixels that 
are equidistant between two or more boundaries. We give preference to the boundary 
sides as follows: horizontal sides are preferred over vertical. The lowest horizontal 
side has the highest priority among horizontal sides, and among all lowest horizontal 
sides, the leftmost has the highest priority. The leftmost vertical side has the highest 
priority among vertical sides, and among all leftmost vertical sides, the highest has 
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the highest priority. In fact, any tie-breaking rule can be applied as long as it is 
applied consistently. The particular choice of the rule only affects the orientation of 
the combs. 

Voronoi Cell Structure. For a pixel p e P let V{p) denote the element of 
dP, closest to p. In case of ties, the tie-breaking rule is applied. Inspired by 
computational-geometry terminology, we call V{p) the Voronoi side of p. We let 
S{p) denote the length of the path from p to the pixel having V{p) as a side (see 
Fig. 7, right). For a boundary side e G dP we let Voronoi(e) denote the (possibly, 
empty) set of pixels, having e is the Voronoi side: Voronoi(e) = {p G P\V{p) = e}. 
We call Voronoi(e) the Voronoi cell of e. The Voronoi cells of the elements of dP 
form a partition of P, which is called the Voronoi decomposition of P. 

A set of pixels C whose dual graph Gc is a straight path or a path with one bend, 
is called a line. Each line C has a root pixel p, which corresponds to one of the two 
leaves of Gc, and a base, e G dP, which is a side of p. 

A (horizontal) comb C is a union of pixels consisting of a set of vertically adjacent 
(horizontal) rows of pixels, with all of the rightmost pixels (or all of the leftmost 
pixels) in a common column. (A vertical comb is defined similarly; however, by our 
tie breaking rules, we need consider only horizontal combs.) A comb is a special type 
of histogram polygon [6]. The common vertical column of rightmost/leftmost pixels 
is called the handle of comb C, and each of the rows is called a tooth. A leftward 
comb has its teeth extending leftwards from the handle; a rightward comb is defined 
similarly. The pixel of a tooth that is furthest from the handle is the tip of the tooth. 
The topmost row is the wisdom tooth of the comb. The root pixel p of the comb is 
either the bottommost or topmost pixel of the handle, and its bottom or top side, 
e G dP, is the base of the comb. See Fig. 1, left. The union of a leftward comb and 
a rightward comb having a common root pixel is called a double-sided comb. 
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Fig. 1: Left: a comb. The base is bold. The pixels in the handle are marked with asterisks, the 
pixels in the wisdom tooth are marked with bullets. Right: Voronoi cells. The sides of dP are 
numbered 1 ... 28 counterclockwise. The pixels in the Voronoi cell of a side are marked with 
the corresponding number. Voronoi cell of side 3 is a comb; Voronoi cells of sides 6, 11, 17, 
25, 28 are empty; cells of sides 1, 7, 10, 18, 24 are lines, comprised of just one pixel; cells of 
the other edges are lines with more than one pixel. 
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An analysis of the structure of the Voronoi partition under our tie breaking rules 
gives: 

Lemma 1. For a side e G dP, the Voronoi cell of e is either a line ( whose dual graph 
is a straight path), or a comb, or a double- sided comb. By our tie-breaking rule, the 
combs may appear only as the Voronoi cells of horizontal edges. The double-sided 
combs may appear only as the Voronoi cells of (horizontal) edges of length 1. 

Let p be a boundary pixel of P, let e G dP be the side of p such that p G 
Voronoi(e). We denote Voronoi(e) by T(p) or T(e), indicating that it is a unique 
tree (a line or a comb) that has p as the root and e as the base. 

Lower Bounds. We exhibit two lower bounds on the cost of an optimal tour, one 
based on the number of pixels and the other on the Voronoi decomposition of the 
domain. At any time let s{R) be the set of pixels of R covered with snow and also, 
abusing notation, the amount of snow on these pixels. Let d{R) = X]pes(_R) ^ip) ■ 

Lemma 2. Let Rbe a subset of P with the snowblower starting from a pixel outside 
R. Then s{R) and d{R) are lower bounds on the cost to clear R. 

Proof. For the snow lower bound, observe that region R cannot be cleared with fewer 
than s{R) snowblower moves because each pixel of s{R) needs to be visited. 

For the distance lower bound, observe that, in order to clear the snow initially 
residing on a pixel p, the snowblower has to make at least S{p) moves. When the 
snow from p is carried to the boundary of P and thrown away, the snow from at most 
D — 1 other pixels can be thrown away simultaneously. Thus, a region R cannot be 
cleared with fewer than d{R) moves. □ 

NP-Completeness. It is known ( [11, 12]) that the Hamiltonian path problem in 
cubic grid graphs is NP-complete. The problem can be straightforwardly reduced to 
SBP. If G is a cubic grid graph, construct an (integral orthohedral) domain P such 
that G = Gp. Since Gp is cubic, each pixel p G P is a boundary pixel, thus, the 
snowblower can throw the snow away from p upon entering it. Hence, SBP on P is 
equivalent to TSP on G, which has optimum less than n + 1 iff G is Hamiltonian 
(where n is the number of nodes in G). The reduction works for any D > 1. 

The algorithms proposed in this paper show that any domain can be cleared using 
a set of moves of cardinality polynomial in the number of pixels in the domain, 
assuming D > 2. Thus, we obtain 

Theorem 1. If D > 2, the SBP is NP-complete, both in the default model and in the 
adjustable throw model, for inputs that are polygonal domains with holes. 

3 Approximation Algorithm for the Default Model 

In this section we give an 8-approximation algorithm for the case when the snow 
can be thrown in all four directions. We first show how to clear a line efficiently 
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with the operation called line-clearing. We then introduce another operation, the 
brush, and show how to clear a comb efficiently with a sequence of line-clearings 
and brushes. Finally, we splice the subtours through each line and comb into a larger 
tour, clearing the entire domain. The algorithm for the default model, developed in 
this section, serves as a basis for the algorithms in the other models. 
Clearing a Line. Let Che a line of pixels; let p and e be its root and the base. We 
are interested in clearing lines for which the base is a boundary side, i.e., e £ dP. 
Let £ = s{C); let the first J pixels of C counting from p be clear. We assume that 
p is already clear (J > 0); the snow from it was thrown away through the side e as 
the snowblower first entered pixel p. Let C\J denote C with the J pixels clear; let 
£ — J = kD + Denote by {C\J)d the first kD pixels of C\J covered with snow; 
denote by the last r pixels on C\J. The idea of decomposing C\J into {C\J)d 
and Cr is that the snow from {C\J)d is thrown away with k "fully-loaded" throws, 
and the snow from Cr is thrown away with (at most one) additional "under-loaded" 
throw. 

We clear line C starting at p by moving all the snow through the base e and 
returning back to p. The basic clearing operation is a back throw. In a back throw 
the snowblower, entering a pixel u from pixel v, throws u's snow backward onto 
V. Starting from p, the snowblower moves along C away from p until either the 
snowblower moves through D pixels covered with snow or the snowblower reaches 
the other end of C; this is called the forward pass. Next, the snowblower makes a 
U-turn and moves back to p, pushing all the snow in front of it and over e; this is 
called the backward pass. A forward and backward pass that clears exactly D units 
of snow is called a D-full pass. 

Lemma 3. For arbitrary D > A the line-clearing cost is at most 2s{C\p)-\-Ad[C\J). 
For D = 2,3 the line-clearing cost is at most 2s{C \ p) + 2d{C\J). If every pass is 
D-full, the cost is Ad{C\J)for D>Aand 2d{C\J)for Z) = 2, 3. 

Proof The clearing cost is c(/: I J) = c{{C\J)d) + c{Cr) = EjLi 2{J-l + iD) + 
2{i - 1) = 2k.J + Dk{k + 1) - 2fc + 2{£ - 1). The snow lower bound of £ \ p is 
s{C\p) ~ L The ^f;ito«ce lowerbound of (£1 J)^! is d{{C\J)D) = j)^i=i{J+ 
i) = kJ + k{kD + l)/2. 
Thus, 

ciC\J)^2.siC\p)+(2 + j^^^^)dmJ)o) 
If every pass is a D-full pass, then c{Cr) = 0. Therefore, c{C\J) = c{{C\J)d) = 

Clearing a Comb. Let C be a comb with the root p, base e, and handle H of length 
H. Let £i . . .£h he the lengths of the teeth of the comb. Since we are interested 

For ease of presentation, we adapt the following convention. For d G {D, [D/2] } and an 
integer w we understand the equahty w = ad + b as follows: b and a are the remainder and 
the quotient, respectively, of w divided by d. 
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in clearing combs for which the base e is a boundary side (e G dP), we assume 
that pixel p is already clear — the snow from it was thrown away through e as the 
snowblower first entered p. 

Our clearing strategy is as follows. While there exists a line C C C rooted at p, 
such that s{C) > D, we perform as many U-fuU passes on C as we can. When no 
such C remains, we call the comb brush-ready and we use another clearing operation, 
the brush, to finish the clearing. 

A brush, essentially, is a "capacitated" depth-first-search. Among the teeth of a 
brush-ready comb that are not fully cleared, let t be the tooth, furthest from the base. 
In a brush, we move the snowblower from p through the handle, turn into t, reach 
its tip, U-turn, come back to the handle (pushing the pile of snow), turn onto the 
handle, move by the handle back towards p until we reach the next not fully cleared 
tooth, turn onto the tooth, and so on. We continue clearing the teeth one-by-one in 
this manner until D units of snow have been moved (or all the snow on the comb has 
been moved). Then we push the snow to p through the handle and across e. This tour 
is called a brush (Fig. 2). 



Fig. 2: Left: a brush-ready comb. The snow is shown in light gray. Center: a brush, D — 4; 
the part of the brush, traveling through the handle, is bold. Right: the comb after the brush. 



Lemma 4. For arbitrary D > A the comb C can be cleared at a cost of at most 
As{C\p)+Ad{C\p) (at most As{C \p) + 2d{C \p) for D = 2,3j. 

Proof. If s{C \ p) < D, then the cost of clearing is just 2s{C \ p), so suppose, 
s{C\p) > D. Let B be the number of brushes used; let B be the set of pixels cleared 
by the brushes. For b = 1 ... B let ti, and t'^ be the first and the last tooth visited 
during the 6th brush. For b G {!...£? — 1} the 6th brush enters at least 2 teeth, so 

tb > t'b — ^b+l- 

Each brush can be decomposed into two parts: the part, traveling through the 
teeth, and the part, traveling through the handle (see Fig. 2). Since each tooth is 
visited during at most 2 brushes, the length of the first part is at most 4 times the size 
of all teeth, that is, 4s (C \ H). The total length of the second part of all brushes is 
2 X]b=i(*b ^ Thus, the cost of the "brushing" is 
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B B 
c{B) < 2^(4 - 1) +4s(C\?^) < 2^i,, + 4s(C\p) -2 (1) 

6=1 6=2 

since ti < H and H > 2 (for otherwise C is a line). 

There are exactly D pixels cleared during each brush 6 £ {0 ... i? — 1}, and each 
of these pixels is at distance at least t^' from the base of the comb. Thus, the distance 
lower bound of the pixels, cleared during brush b, is at least t^' ■ Consequently, the 
distance lower bound of B 

B B-1 B 

d{B) >^h,>Y^ h+i =^tb (2) 

6=1 6=1 6=2 

From (1) and (2), B can be cleared at a cost of at most 2d{B) + 4s(C \ p). 

Let V Q Che the pixels, cleared during the line-clearings. By our strategy, during 
each line-clearing, every pass is Z?-full; thus, by Lemma 3, V can be cleared at a cost 
of at most 4(i(7') (or 2d{V) if D = 2, 3). Since V and B are snow-disjoint and 
V yj B = C\p,the, lemma follows. □ 

The above analysis is also valid in the case when the handle is initially clear. This 
is the case when the second side of a double-sided comb is being cleared. Thus, a 
double-sided comb can be cleared within the same bounds on the cost of clearing. 

Clearing the Domain. Now that we have defined the operations which allow us to 
clear efficiently lines and combs, we are ready to present the algorithm for clearing 
the domain. 

Theorem 2. For arbitrary D > A (resp., D = 2,3j an 8-approximate (resp., 6- 
approximate) tour can be found in polynomial time. 

Proof. Let pi, . . . , pM be the boundary pixels of P as they are encountered when 
going around the boundary of P counterclockwise starting from g = pi; let 
ei, . . . , Cm G dP be the boundary sides of pi, . . . ,p]^f such that e,; = Ve{pi), i = 
1 . . . M. The polygon P can be decomposed into disjoint trees T{pi), . . . , T{pm) = 
T(ei), . . . , T{eM) with the bases ei . . . , cm, where each tree T(ei) is either a line 
or a comb. 

Our algorithm clears P tree-by-tree starting with T(ei) = T{g). By Lem- 
mas 3 and 4, for i = 1 . . . M, the tree T{pi) \ pi can be cleared at a cost of at 
most 4s(T(pi) \ Pi) + Ad{T{pi) \ pi) starting from pi and returning to pi. Since 
T{pi) \ Pi = P \ {pi . . .pAi}, the interior of P can be cleared at a cost 
of at most c{P \ {pi . . .pm}) ^ 4s(P \ {pi . . .pAi}) + Ad{P \ {pi . ..pAi}) < 
4s(P \g) + 4:d{P \ g) — AM + 4. Finally, the tours clearing the interior of P can be 
spliced into a tour, clearing P at a cost of at most 2M. Since the optimum is at least 
s{P \ g) and is at least d{P \ g), the theorem follows. □ 

4 Other Models 

In this section we give approximation algorithms for the case when the throw 
direction is restricted. Specifically, we first consider the adjustable throw direction 
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formulation. This is a convenient case for the snowblower operator who does not 
want the snow thrown in his face. We then consider the fixed throw direction model, 
which assumes that the snow is always thrown to the right. 

We remark that the low approximation factors of the algorithms for the default 
model, presented in the previous section, were due to a very conservative clearing: 
the snow from every pixel p € P was thrown through the Voronoi side V{p). Un- 
fortunately, it seems hard to preserve this appealing property if throwing back is 
forbidden. The reason is that the comb in the Voronoi cell Voronoi(e) of a boundary 
side e £ dP often has a "staircase"-shaped boundary (see Fig. 8); clearing the first 
"stair" in the staircase cannot be done without throwing the snow onto a pixel of 
Voronoi(e'), where e' 7^ e is another boundary side. That's why the approximation 
factors of the algorithms in this section are higher than those in the previous one. 

Adjustable Throw Direction 

In the adjustable-throw model the snow cannot be thrown backward but can be 
thrown in the three other directions. To give a constant- factor approximation algo- 
rithm for this case, we show how to emulate line-clearings and brushes avoiding back 
throws (Fig. 3). The approximation ratios increase slightly in comparison with the 
default model. 

Line-clearing. We can emulate a (half of a) pass by a sequence of moves, each with 
throwing the snow to the left, forward or to the right (Fig. 3, left and center). Thus, 
the line-clearing may be executed in the same way as it was done if the back throws 
were allowed. The only difference is that now the snow is moved to the base when 
the snow from only [15/2 J pixels (as opposed to D pixels) of the line is gathered. 

Lemma 5. The line-clearing cost is at most 3D/ \_D/2\d{C\J) + 2s(£ \p). If every 
pass is [D/2\-full, the cost is iD /\_D /2\d{C\J). 

Proof. Let £ - J = fc' [D/2\ + r' . Let {C\J)yd/2\ be the first k'\_D/2\ pixels of 
£| J, let Cr' be its last r' pixels. Then the cost of the clearing of C\J is c{C\J) = 

c{{C\J)yD/2\)+c{Cr:) = J2ti2{J+i[D/2\)+2£ ^ 2k' J+lD /2\k'{k' + l) + 2£. 
The lower bounds are given by s{£\p) = £ — 1 and 




i=l 



(3) 



Thus, 



c{C\J) < 




2 + 



k'J+^k'^ + !^ 



d{C\p) + 2s{C\p) 



□ 



Brush. Brush also does not change too much from the default case. The difference 
is the same as with the line-clearing: now, instead of clearing D pixels with a brush. 
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Fig. 3: Emulating line-clearing and brush. The (possible) snow locations are in light gray; s is 
the snowblower. Left: forward and backward passes in the default model; there are D units of 
snow on the checked pixel. Center: the passes emulation; there is (at most) 2[_D/2J units of 
snow on the checked pixel. Right: the snow to be cleared during a brush is in light gray; there 
is [-D/2J light gray pixels. Together with each gray pixel, the snow from at most 1 other pixel 
is moved — thus, the brush is feasible. 



only [D/2\ pixels are (considered) cleared (Fig. 3, right). The definition of a brush- 
ready comb is changed too — now we require that there is less than [15/2 J pixels 
covered with snow on each tooth of such a comb. 

Lemma 6. A comb can be cleared at a cost ofiD /^D /2\d{C \ p) + 4s(C \ p). 

Proof. In comparison with the default model (Lemma 4) several observations are in 
place. The number of brushes may go up; we still denote it by B. The cost of the 
brushes 1 ... B — 1 does not change. If the Bth brush has to enter the first tooth, 
there may be 2 more moves needed to return to the root of the comb (see Fig. 3, 
right); hence, the total cost of the brushing (1) may go up by 2. The distance lower 
bound (2) goes down by D/ [D/2J . The rest of the proof is identical to the proof of 
Lemma 4 (with Lemma 5 used in place of Lemma 3). □ 

Observe that in fact the snow can be removed from more than \D/2\ pixels during 
a brush; we just ignore it for now in our analysis. Note that a double-sided comb can 
also be cleared in the described way. 

Clearing the Domain. As in the default case (Theorem 2), 

Theorem 2>. A (4 + 3D / \D /2J )-approximate tour can be found in polynomial time. 

Comment on the Parity of D. We remark that if D is even, the cost of the clearing is 
the same as it would be if the snowblower were able to move through snow of depth 
D + 1 (the slight increase of 6/(1? — 1) in the approximation factor would be due to 
the decrease of the distance lower bound). 

Fixed Throw Direction 

In reality, changing the throw direction requires some effort. In particular, a snow 
plower does not change the direction of snow displacement at all. In this section we 
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consider the fixed throw direction model, i.e., the case of the snowblower which 
can only throw the snow to the right. We exploit the same idea as in the previous 
subsection — reducing the problem in the fixed throw direction model to the problem 
in the default model. All we need is to show how to emulate line-clearing and brush. 
In what follows we retain the notation from the previous section. 

Lemma 7. The line-clearing cost is at most 24D/[Z)/2jd(£| J) + 25s(£ \ p). If 
every pass is [D/2\-fiill, the cost is 2AD /[D l2\d{C\J). 

Proof. We first consider clearing a line whose dual graph is embedded as a single 
straight line segment and whose base is perpendicular to the segment; we describe 
the line-clearing, assuming that the line is vertical. Next, we extend the solution to 
the case when the base is parallel to the edges of the dual graph; this can only be a 
horizontal line — the first tooth in a (double-)comb. Finally, we consider clearing an 
L-shaped line; this can only by a tooth together with the (part of the) handle. 

A Line C with Gc^e. As in the adjustable-throw case (see Fig. 3, left and center), 
to clear C we will need to use the pixels to the right of C to throw the snow onto. 
Let p' be the pixel, following p in dP. Before the line-clearing is begun, it will be 
convenient to have p' clear. Thus, the first thing we do upon entering C (through p) is 
clearing p'. Together with returning the snowblower to p it takes 2 or 4 moves (Fig. 4, 
left); we call these moves the double-base setup. 



m 



Fig. 4: Left: the double-base setup. Right: before the forward pass the snow below the snow- 
blower is cleared on both lines. 



Then, the following invariant is maintained during line-clearing. If the snow- 
blower is at a pixel q £ C before starting the forward pass, all pixels on C from p 
to q are clear, along with the pixels to the right of them (Fig. 4, right). The invariant 
holds in the beginning of the line-clearing and our line-clearing strategy respects it. 

Each back throw is emulated with 5 moves (Fig. 5, left). After moving up by 
[-D/2J pixels (and thus, gathering 2[_D/2J units of snow on these [-D/2J pixels), the 
snowblower U-turns and moves towards p "pushing" the snow in front of it; a push 
is emulated with 1 1 moves (Fig. 6). 

The above observations already show that the cost of line-clearing increases only 
by a multiplicative constant in comparison with the adjustable-throw case. A more 
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careful look at the Figs. 5, 6 reveals that: (1) in the push emulation the first two moves 
are the opposites of the last two, thus, all 4 moves may be omitted - consequently, 
a push may be emulated by a sequence of only 7 moves; (2) if the boundary side, 
following e, is vertical, the last push, throwing the snow away from P, may require 9 
moves (Fig. 5, right); and, (3) when emulating the last back throw in a forward pass, 
the last 2 of the 5 moves (the move up and the move to the right in Fig. 5, left) can 
be omitted - indeed, during the push emulation, the snowblower may as well start 
to the right of the snow (see Fig. 6). Thus, a line C\J can be cleared at a cost of 

cmj) < 4 +Eti - 1 + - 1) L#J + 5 LtJ + nj + LtJ - 1)) + ^+ 

5r' + 7{£ - 1). 



Fig. 5: Left: emulating back throw. Right: pushing the 2 [-D/2J units of snow away from P 
and returning the snowblower to p may require 9 moves. 



■ 
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Fig. 6: Emulating pushing the snow in front of the snowblower. 



A Line C with G£||e. Consider a horizontal line, extending to the left of the base; 
such a line may represent the first tooth of a comb. The double-base can be cleared 
with 8 or 12 moves (see Fig. 9), the root can be cleared with 3 moves (see Fig. 10, 
left) instead of 9 moves (see Fig. 5, right); the rest of the clearing does not change. 

Consider now a horizontal line extending to the right of the base (see Fig. 10, 
right); such a line may appear as the first tooth in a double-sided comb. The double- 
base for such a line can be cleared with 3 moves (see Fig. 10, right); the rest of the 
clearing is the same as for the vertical line. 

L-shaped Line. An L-shaped line C consists of a vertical and a horizontal segment. 
Each of the segments can be cleared as described above. 

Thus, any line C\J can be cleared at a cost of at most c{C\J) < 12 + 

Eti ( J - 1 + - 1) Lf J + 5 Lf J + 7( J + z Lf J - 1)) + J + 5r' + 7(£ - 1). 
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Since the snow and distance (3) lower bounds do not change, the lemma follows. 

□ 

Lemma 8. A comb can be cleared at a cost o/34s(C \ p) + 24D/ [_D/2J d(C \ p). 

Proof. Brush in the fixed throw direction model can be described easily using anal- 
ogy with: a) brush in default and adjustable-throw models and b) line-clearing in 
fixed-throw model. As in the adjustable-throw model, we prepare to clear [-D/2J 
pixels during each brush. Same as with line-clearing, we setup the double-base for 
the comb with at most 12 moves; also, 9 moves per brush may be needed to push 
the snow away from P through the base. Back throw and push can be emulated with 
5 and 7 moves (Fig. 5, left and Fig. 6). Thus, if the cost of a brush (1) in the de- 
fault model was, say, c, the cost of the brush in the fixed-throw model is at most 
7c + 9. Since any brush starts with the double-base setup, c > 6; this, in turn, implies 
7c+9 < (51/6)c. Hence, the cost of clearing B increases by at most a factor of 51/6. 

By Lemma 7, the cost of clearing V, c{V) < 2AD /[D /2\d{V). The snow and 
distance lower bounds do not change in comparison with the adjustable-throw case. 
The lemma now follows from simple arithmetic. □ 

As in the default and adjustable-throw models (Theorems 2,3), 

Theorem 4. A (34 + 2AD / [_D/2J ) -approximate tour can be found in polynomial time. 

5 Extensions 

Polygons with Holes. Our methods extend to the case in which P is a polygonal 
domain with holes. There are two natural ways that holes may arise in the model. 

In the first, the holes represent obstacles for the snowblower (e.g., walls of build- 
ings that border the driveway). No snow can be thrown onto such holes; the holes' 
boundaries serve as walls for the motion of the snowblower and for the deposition of 
snow. Our algorithm for the default model extends immediately to this variation. The 
SBP in restricted-throw models, however, may become infeasible; it is not always 
possible for the snowblower to enter a thin (one-pixel-width) "channel" between two 
walls. 

In the second variation, the holes' boundaries are assumed to be the same "cliffs" 
as the polygon's outer boundary. It is in fact this version of the problem that we 
proved to be NP-complete. With some modifications our algorithms work for this 
variation as well; we sketch here the necessary modifications. The boundary, dP, 
of the domain now consists of the boundary pixel sides both on the outer boundary 
of the polygon, denoted d°P, and on the boundaries of the holes. As in the simple 
polygon case, we build the Voronoi decomposition of the domain and prepare to 
clear it Voronoi-cell-by-Voronoi-cell. The only problem with it now is that there is no 
readily available Hamiltonian cycle through the elements of dP. To build a suitable 
cycle we find a spanning tree, T, in the graph having a node for each hole and a node 
for d"P and having an edge between two holes Hi and H2 if there are boundary 
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pixel sides, ei e dHi and e2 G 9i?2, such that there exist adjacent pixels pi G 
Voronoi(ei) andp2 G Voronoi(e2). Then, after Voronoi(ei) is cleared, we direct the 
snowblower to start clearing Voronoi(e2).^ After that, the snowblower continues to 
clear the Voronoi cells of the boundary pixel sides of H2 until another edge of T, 
connecting H2 to a hole H3, is encountered (if H2 is a leaf of T, H3 = Hi), and so 
on. The cost of clearing the Voronoi cells of the sides in dP does not change, and 
the cost of traversing T is at most twice the size of P; thus, the proposed algorithm 
remains a constant-factor approximation algorithm. 

Nonrectilinear Polygonal Domains. Our discussion so far has assumed that P is 
integral-orthogonal; in this case, the snowblower makes axis-parallel movements. If 
P is rectilinear, but not necessarily integral, we proceed as in [3]: first, the bound- 
ary of P is traversed once, and then our algorithms are applied to the remaining 
part, P', of the domain P. Every time the snow is thrown away from P', a certain 
length (which depends on the throw model) may need to be added to the cost of the 
tour; thus, the approximation factors of our algorithms may increase by an additive 
constant. 

We can also extend our methods to general nonrectilinear domains. Since the 
snowblower is not allowed to move outside the domain, care must be taken about 
specifying which portion of the domain is actually clearable. This portion can be 
found by traversing the boundary of the domain; then, the accessible portion can be 
cleared as described above. 

Implicit Representation of the Tour. As in [3,4], we make the distinction be- 
tween explicit and implicit representations of the domains and snowblower tours. As 
mentioned in the introduction, we assume that a domain is given as the union of pix- 
els; this way the size of the input to the problem is 0(N), where N is the number 
of pixels in the domain. The size of the description of the snowblower tour produced 
by our algorithm is polynomial in N, i.e., polynomial in the input size. Instead, the 
domain may be given in polygonal representation, as a list of coordinates of its n 
vertices; the size of such a representation is 0{n\ogW), where W is the largest 
coordinate in the input. In principle, TV may be n{W), and hence the length of the 
description of the snowblower tour may appear to be exponential in the size of the 
input. To avoid such a blow-up in the size of the tour, it is possible to represent the 
tour succinctly so that the size of the representation is polynomial in n log W. 

Indeed, our algorithms produce tours, comprised of line-clearings and brushes. 
Our Voronoi decomposition of the domain is the discretized version of the Voronoi 
diagram of the edges of the domain; the latter is 0{n) in size and can be found 
efficiently ( [5]). Given the diagram, it is easy to constrain the (axis-parallel) motion 
of the snowblower to stay within a Voronoi cell of an edge: when "in doubt", i.e., 
when the snowblower is about to enter a pixel, intersected by an edge of the Voronoi 
diagram, it can be decided "in place", in constant time (based on the tie breaking 
rules), whether entering the pixel will place the snowblower into the Voronoi cell of 
another side of P. 

^ The clearing starts from creating a path of width 2 from p2 to 62; the path "bridges" the 
holes Hi and H2- 
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The clearing of a Voronoi face in the Voronoi diagram is done tree-by-tree: first 
a (double-)comb is cleared (if present), then a set of lines (if present), then the other 
comb (if present). Thus, finding a short representation of a tour boils down to exhibit- 
ing succinct representations of the tours through a comb and through a set of lines 
with adjacent bases comprising a boundary edge of P. The descriptions of these tours 
given in Sections 3 and 4 provide such representations. 

The Vacuum Cleaner Problem. Consider the following problem. The floor - a 
polygonal domain, possibly with holes - is covered with dust and debris. The house 
is equipped with a central vacuum system, and certain places on the boundary of 
the floor (the baseboard) are connected to the "dustpan vac" - a dust dump location 
of infinite capacity [1]. The robotic vacuum cleaner has a dust/debris capacity D 
and must be emptied to a dump location whenever full. The described problem is 
equivalent to the SBP in default throw model, and actually provided the motivation 
to study the SBP with throwing backwards allowed. 

Nonuniform Depth of Snow. Our algorithms generalize straightforwardly to the 
case in which some pixels of the domain initially contain more than one unit of 
snow. For a problem instance to be feasible it is required that there is less than D 
(less than [I3/2J in restricted-throw direction models) units of snow on each pixel. 
The approximation ratios in this case depend (linearly) on D (or, in general, on the 
ratio of D to the minimum initial depth of snow on P). 

Capacitated Disposal Region. Another generalization of the problem is to con- 
sider the dump locations to have finite capacities. If instead of "cliffs" at the bound- 
ary of P, there is a finite capacity (maximum depth) associated with each point in 
the complement of P, the SBP more accurately models some material handling prob- 
lems, but also becomes considerably more difficult. The snow lower bound still ap- 
plies, the distance lower bound transforms to a lower bound based on a minimum- 
cost matching between the pixels in P and the pixels in the complement of P. This 
problem represents a computational problem related to "earth-mover distance" [7] 
and is beyond the scope of this paper 

Open Problems. The complexity of the SBP in simple polygons and the complex- 
ity of the SBP in the fixed-throw model are open. We also do not have an algorithm 
for the case of holes-obstacles and restricted throw direction; the hardness of this 
version is also open. 
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Appendix: Figures 




Fig. 7: Left: the tie breaking tree; on a tie, the upper edge of the tree is followed. Right: an 
example; by the tie breaking rules, V{p) — e, the bold edge; 5{p) = 4. 
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Fig. 8: Part of the comb of Voronoi(S). Ve{p) C E. Ve{a),Ve{b),Ve{c) ^ E, so the 
snow form a must be thrown onto Voronoi(_E') 7^ Voronoi(£'). Nevertheless, observe that 
S{a) = S{p) — 1 - one could try to make use of it. 
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Fig. 9: Setting up the double-base for clearing a horizontal line extending to the left of its 
base. Depending on the direction of the edge adjacent to the base from the right, there are 8 
(above) or 12 (below) moves necessary. 
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Fig. 10: Left: clearing the root of a horizontal, extending to the left, line with 3 moves. There 
is 2 [-D/2J units of snow on the checked pixel. Right: Setting up the double-base for clearing 
a horizontal line extending to the right of its base. 



